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ABSTRACT 

During  the  design  stage  of  automatic  control  systems  for 
aircraft,   fire  control  or  other  physical  systems,   system  com- 
ponents must  be  selected.     Values  of  the  parameters  of  these 
components  must  be  determined  so  as  to  produce  the  desired 
dynamic  properties  of  the  resultant  system. 

The  method  of  estimating  dynamic  characteristics  employ- 
ed in  this  investigation  is  based  on  the  encirclement  criteria  of 
Cauchy's  Residue  Theorem,.       A  computer  was  developed  to  aid  in 
this  estimation  by  making  certain  modifications  to  the  data  pre- 
sentation of  the  Fourier  Harmonic  Synthesizer. 

This  investigation  indicated  that  dynamic  characteristics 
may  be  readily  estimated  by  the  method  to  be  described  for  any 
selected  set  of  system  parameters.     System  parameter  specifi- 
cations and  limits  may  also  be  set  when  the  dynamic  character- 
istics are  specified. 

Several  examples  are  included  to  illustrate  the  techniques 
involved  and  to  indicate  the  suitability  of  applying  this  method  to 
aid  in  the  analysis  of  higher  order  systems. 
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Title:  Assistant    Professor 

Aeronautical  Engineering 
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OBJECT 


The  object  of  this  investigation  is  to  present  suitable 
techniques,   based  on  the  Cauchy  Residue  Theorem  and  utilizing 
a  modified  Fourier  Harmonic  Synthesizer  as  a  computer,  to  aid 
in  the  estimation  of  the  dynamic  characteristics  of  physical 
systems. 
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CHAPTER  I 
INTRODUCTION 


In  the  design  of  physical  systems  some  of  the  most  important 
considerations  are  the  dynamic  properties  of  the  resultant  system. 
It  will  be  shown  in  Summary  1  that  the  differential  equation  des- 
cribing the  system  may  be  written  with  the  input  terms  on  the 
right  side.     By  setting  the  right  side  equal  to  zero  the  homogeneous 
equation  is  obtained.     This  equation  determines  thetransient  modes 
of  the  system,  which  are  independent  of  the  input.     The  actual  in- 
put affects  the  magnitude  of  the  transient  modes  but  has  no  effect 
on  the  damping  ratios,  the  oscillation  frequencies,   or  response 
times  of  the  transients.     These  characteristics  are  determined 
by  the  component  parameters.     Once  the  system  is  built,   control 
over  these  dynamic  properties  is  restricted.     Since  these  proper- 
ties have  a  prominent  effect  on  the  performance  of  the  system,  the 
engineer  must  have  some  means  of  predicting  them  during  the 
design  stage. 

Exact  determination  of  system  parameters  is  often  neither 
necessary  nor  desirable  in  the  early  stages  of  design  if  limiting 
values  of  the  parameters  can  be  estimated  in  terms  of  the  dynamic 
characteristics  of  the  system.     Such  an  estimate  will  enable  the 
engineer  to  set  limits  on  the  component  specifications  in  the  in- 
itial stages  of  design. 

The  method  of  estimation  of  system  dynamic  characteristics 
described  in  this  thesis  is  based  on  a  concept  first  advanced  by 
Sidney  Lees  in  Massachusetts  Institute  of  Technology  Instrumen- 
tation Laboratory  Report  R-71  of  April  1954   *  '  . 

The  Mapping  Theorem,  described  in  Instrument  Engineer- 
ing Volume  II  '    '  and  in  Theory  of  ServoMechanisms,   V3)    and 
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the  Cauchy  Residue  Theorem  described  by  Guillemin  *   '   and 
others,   are  the  basis  for  this  concept. 

Lees  applies  this  theorem  '  '  and  describes  a  procedure  by 
which  dynamic  characteristics  may  be  estimated  from  polynomial 
contours  approximated  by  the  computation  of  certain  selected 
points. 

If  these  contours  can  be  accurately  computed  and  displayed 
by  some  device,  the  application  of  the  Cauchy  Theorem  to  this 
problem  should  be  greatly  facilitated. 

A  Fourier  Harmonic  Synthesizer  built   at  the  Massachusetts 
Institute  of  Technology  Instrumentation  Laboratory  and  described 
by  Instrumentation  Laboratory  Engineering  Memo  6445-E-42  '    ' 
was  used  for  this  purpose.  Modifications  to  the  existing  instal- 
lation and  the  techniques  employed  to  utilize  this  device  as  an 
aid  in  the  estimation  of  the  dynamic  characteristics  of  physical 
systems  are  described  in  the  succeeding  chapters. 
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CHAPTER  2 
MATHEMATICAL  BACKGROUND 


The  performance  of  physical  systems  may  often  be  described 
by  linear  differential  equations  with  coefficients  that  are  constant 
once  they  have  been  decided  upon.     The  properties  of  the  compo- 
nents of  the  system  may  be  obtained  from  the  equation  coefficients. 

The  complete  solution  of  the  linear  differential  equation 
consists  of  the  transient  solution,  plus  the  forced  solution  as  is 
shown  in  Summary  1. 

The  transient  solution  of  the  linear  differential  equation  is 
the  solution  of  the  homogeneous  equation,  which  is  obtained  by 
setting  the  right  hand  side  of  the  differential  equation  equal  to  zero. 

SUMMARY  I 

As  expressed  in  Instrument  Engineering  Volume  II  !  '    '  ,   a 
typical  ordinary  linear  differential  equation  has  the  form 

n  m  . 

dlV  V     u         du  n. 

ai      ~T       =       >      bk    k~  (l) 

1     dt1  L       K     dtk 

i=0  k=0 

The  homogeneous  equation  is 

i  =  0 

The  associated  characteristic  equation  is 
n 

Y     a.     z1    =    0  (3) 

i=0 
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A  fundamental  theorem  of  algebra  described  in  Summary 
24-1  Instrument  Engineering  Vol.     II  (2)  states  that  every 
equation  of  nth  degree  has  n  and  only  n  roots.     The  roots  of 
Equation  (3)  are  z,,   z„,   z„,    .  .  .  .  z    „     In  general  the  roots  are 
complex  numbers  so  that 

zi  =  xi  +  iYi  (4) 

When  all  roots  are  distinct,  transient  solutions  of  Equation 
(1)  have  the  form 

n  ,  n  ,,  , 

z.t  x.t       jy.t 

V(tr)    =  >       ci   e  =     I     ci   e       l    e  '  (5) 

i=l  i=l 

Roots  of  the  characteristic  equation  must  have  the  di- 
mensions of  inverse  time  (from  Equation  (5)  ).     Real  roots  may 
be  expressed  as 

Zi  =  Xi=  -(CT}.  =  if-  <6> 

where  (CT).  =  t.  =  characteristic  time  of  real  root. 
Complex  roots  may  have  the  form 


zk  =  xk  +  ^k  =  "  *k  wn,    +  J   wn,    J~l  ~  Ik  <7> 


where  £,,=  (DR),    =  damping  ratio  of  kth  root, 


u       =  undamped  natural  frequency  of  kth  root, 
k 
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_L_    =  ._L       —  »■ ± 

*k  *JT)k  Tk 

characteristic  time  of  kth  root. 


andxk;     ■  ^k  V  (CT),        ~     ^ere  (CT)k         Tl 


In  polar  coordinates,   the  kth  complex  root  form  is 
zk  =  xk  +  ^k  =   V    e   J0=un,    e    J   [*"(ADR)k] 


J  (*  -  A      ) 
u        e  ^k  (8) 


nk 


where      9   =    ?r    -  (ADR),     =    vr   -   A 


k 


1 


and  (ADR),    =  A„     =  cos      £,     =  angle  of  damping  ratio  of  the 
kth  root. 

The  transient  solution  determines  what  oscillatory  modes 
can  be  present,   and  their  mode  shapes,   and  also  determines  the 
extent  to  which  the  dynamic  performance  of  the  system  will  be 
affected  by  the  specified  input. 

The  forced  solution  reflects  the  effect  of  the  input  and 
determines  the  relative  magnitude  of  the  modes  indicated  in  the 
transient  solution  superimposed  on  the  forced  solution.     It  also 
determines  the  system  response  time  due  to  the  effect  of  the  in- 
put on  the  mode  shape  but  has  no  effect  on  the  damping  or  oscil- 
lation frequencies  of  the  second  order  modes.     The  response  time 
for  the  system  will,   in  no  case,   be  greater  than  that  indicated  by 
the  slowest  decaying  term  of  the  transient  solution. 

This  investigation  is  not  restricted  to  any  specified  input 
and  considers  only  dynamic  effects  resulting  from  the  roots  of  the 
characteristic  equation.     The  analysis  of  this  part  of  the  complete 
performance  equation  gives  information  on  dynamic  properties 
which  greatly  affect  the  performance  of  the  system  for  any 
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arbitrary  input  and  it  is  to  this  end  that  the  investigation  is 
directed. 

The  roots  of  the  homogeneous  equation  expressed  as  a 
characteristic  equation  with  the  variable  z  are  defined  in  M,  I.  T. 
Instrumentation  Laboratory  Report  R-71  0)  as  the  dynamic 
characteristics  of  the  associated  system.     The  dynamic  charac- 
teristics can  be  determined  if  the  roots  of  the  characteristic 
equation  can  be  found  by  factoring  or  other  means. 

Lees  (^'discusses  the  important  considerations  in  the 
investigation  of  transient  solutions  of  most  practical  systems. 
One  consideration  is  that  the  transient  solutions  decay  to  zero, 
thus  providing  absolute  stability.     In  addition,   it  is  necessary 
that  the.  transient  solutions  decay  rapidly  enough  to  meet  specified 
conditions;   this  is  the  relative  transient  stability  requirement. 
The   damping  ratio  and  the  oscillation  frequency  are  also  impor- 
tant considerations  in  investigating  the  transient  solutions. 
These  are  the  dynamic  characteristics  and  are  apparent  from  the 
transient  solution  or  roots  of  the  characteristic  equation.       In  this 
thesis,  the  terms   t,    £,  and  to     specified  by  the  roots  of  the  charac- 
teristic equation  a*re  considered  to  be  the  dynamic  characteristics 
of  the  physical  system.     The  equivalence  of  these  terms  and  the 
characteristic  equation  roots  are  shown  in  Figure  1. 


x-     l-      l 


Real  Root 


Imaginary 
Axis    (y) 


ADR=A 


Real  Axis  (x) 


Complex  Roots 


Figure  1:    Roots  of  the  Characteristic  Equation  as  Represented 

on  the  CornpLex   Plane.. 
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Absolute  transient  stability  of  the  system  exists  when  all 

of  the  roots  of  the  characteristic  equation  have  negative  real  parts, 

x.t  x,  t 

*  IK. 

i.e.  ,  when  x.  and  x,    are  negative,   e  and  e  decay  to  zero  as 

1  XV 

t  —  oo. 

Relative  transient  stability  exists  when  the  smallest  of  the 
magnitudes  of  the  real  parts  of  the  roots  is  greater  than  a 
specified  reference  magnitude,   i.  e.  ,   if 


x.    * 

1 
mp 


x      J  then 
ref 


(CT).   <  (CT)      .    (x.  is  most  positive  real  part  of  all  the  roots). 

mp 

This  has  been  discussed  in  greater  detail  in  Instrument  Engineering, 
Volume  II  (2\  Chapter  20. 

The  damping  ratio  and  oscillation  frequency  of  the  system 
are  related  to  the  real  and  imaginary  values  of  the  complex  roots 
as  shown  in  Equation  (7). 

In  many  practical  systems  the  dynamic  characteristics  are 
restricted  to  a  certain  desired  sector  or  range  depending  on  the 
relative  transient  stability  requirements  and  specified  undamped 
natural  frequency  and  damping  ratio  limitations.     These  sectors 
and  ranges  are  shown  in  Figure  2  below.     Specification  of  a  sector 
imposes  limitations  on  the  roots  and  thus  the  coefficients  Of  the 
characteristic  equation,   and  it  is  desired  to  determine  the  range 
of  each  coefficient  which  will  provide  dynamic  characteristics 
within  specified  limits.     Figure  2  shows  how  these  limitations  may 
be  investigated  with  respect  to  the  characteristic  equation  complex 
plane  root  loci.     These  contours  are  useful  when  investigating  for 
relative  stability,  oscillation  frequency  and  damping  ratio.     In 
most  practical  systems  the  specifications  impose  limitations  on 
two  or  all  three  of  these  characteristics  and  a  contour  similar  to 
Figure  2-d  is  used. 
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£Wj  j 


(a) 

(b) 

(c) 

(d) 

Relative 

Undamped 

Damping 

Superimposing  Figures 

Transient 

Natural 

Ratio 

a.    b.    and  c 

Stability 

Frequency 

Limits 

Dynamic  C'harac  teristii 

Limit 

Limits 

Limits 

Figure  2:    Standard  Contours 

This  investigation  is  concerned  with  the  techniques  con- 
cerning the  solution  of  two  problems:   determination  or  estimation 
of  the  dynamic  characteristics,   and  fixing  the  limitations  and  re- 
lationships of  the  coefficients.     Fixation  of  the  coefficients  of  the 
characteristic  equation  determines  the  limits  of  the  component 
parameters  in  the  actual  system. 

The  simpler  case,  that  of  locating  the  roots  of  the  charac- 
teristic equation  with  known  constant  coefficients  will  be  taken  up 
first.     This  case  will  be  treated  in  some  detail  since  all  the 
principles  involved  are  equally  applicable  in  estimating  the  dy- 
namic characteristics  from  equations  with  adjustable  coefficients. 

In  this  presentation,  the  methods  used  in  estimating  the  dy- 
namic characteristics  of  physical  systems  are  based  on  Cauchy's 
Residue  Theorem,   as  presented  by  Guillemin  (**)  and  others. 

From  Equation  (3), 


n 


a.    z 
1 


i  =  0 


20 


This  is  expanded  as 


P(z)  =  a     zn+a       ,  z  n  l+  .  .  .+  a„  z     +  a,  z  +  a      =0  (9) 

v/n  n-1  2  I  o  v   ' 


From  Cauchy's  Residue  Theorem,   if  P(z)  is  a  polynomial 
of  z,    P'  (z)  is  the  derivative  of  the  polynomial  with  respect  to 
z,   and  C  is  a  closed  countour  on  the  z  complex  plane  containing 
N  roots  of  P(z)  =  0, 


r 

i 


¥$■  dz  =  j  2tn 


(10) 


This  integral  may  be  interpreted  to  mean  that  the  argument 
of  P(z)  increases  exactly  2tfN  times  as  the  closed  contour  C  is 
traced  out  in  a  positive  sense,  or  that  the  P(z)-plane  path  encircles 
it  u  origin  N  times.     This  is  also  referred  to  as  the  mapping  theorem, 
in  Instrument  Engineering,   Volume  II  (2),   and  is  illustrated  in 
Figure  3. 


Closed  contour 
C  in  z- plane 


(a)    Three  roots  on  z-plane 
within   contour 


P(z)-plane  path 
corresponding 
to  C 


(b)    Three  encirclements  of  the 
origin  on  the  P(z)=plane 


Figure  3.    Illustration  of  Mapping  Theorem 
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A  closed  contour  similar  to  Figure  2-d  and  corresponding 
to  Figure  3~a  may  be  drawn  on  the  characteristic  equation  root 
complex  plane,    in  which  it  is  desired  to  investigate  for  roots  of 
the  characteristic  equation.        By  substituting  values  of  the 
variable  lying  along  the  complex  plane  contour  of  Figure  2~d,   a 
corresponding  contour  on  the  P  (z)  plane  corresponding  to  Figure 
3-b  may  be  determined.     The  number  of  encirclements  of  the 
P(z)  plane  origin  is  the  number  of  polynomial  roots  within  the 
root  complex  plane  contour.     As  will  be  shown  in  the  examples, 
the  polynomial  P(z)  of  Equation  (9)  usually  includes  a  constant 
term  a   .     It  is  convenient  to  write  the  polynomial  as  P(z)-a 
and  -a     becomes  the  P(z)  plane  encirclement  point. 

The  P(z)  plane  path  of  Figure  3-b  may  be  calculated  by 
hand  by  substituting  values  of  the  variable  along  the  contour  of 
Figure  2-d.     This  becomes  impractical  since  a  large  number  of 
calculations  is  required  if  a  complete  and  accurate  P(z)  plane 
path  corresponding  to  the  z  plane  contour  of  Figure   2=d  is  desired. 
Chapter  3  describes  how  the  Fourier  Synthesizer  may  be  used  to 
obtain  the  information  provided  by  the  P(  z)  plane  path. 
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CHAPTER  3 


APPLICATIONS  OF  THE  FOURIER  HARMONIC  SYNTHESIZER 


The  Fourier  Harmonic  Synthesizer  may  be  used  to  facili- 
tate the  application  of  Cauchy's  Residue  Theorem  to  this  problem, 
reduce  the  required  calculations,   increase  the  accuracy  and  bring 
other  related  items  of  interest  to  light. 

It  has  been  shown  that  the  polynomial  in  z   may  be  written 

in  the  form 

■d/v  ^n         in0,  T-,n-l         i(n-l)    5  , 

P  (z)  -  a     =  a     R       eJ        +  a       ,    R  eJV         '        +   .  .  . 

v   '        o        n  n  -I 


+  a ,      Re 


jo  m 


by  setting  z  =  Re  ■"    . 

This  polar  form  of  the  equation  is  suitable  for  setting  into 
the  harmonic  synthesizer.     The  method  of  scaling  to  obtain  full 
scale  deflection  of  the  synthesizer  and  maximum  accuracy  there- 
by is  given  and  explained  in  Example  I. 

Each  of  the  coefficients   (R      ,   R    "    )    corresponds  to  the 
radius  of  the  contribution  of  the  particular  term  of  each  order  as 
the  variation  is  traced  out  from  6   =  0  to  9  -  2  it. 

The  synthesizer  traces  the  polynomial   I  P  (zuplane    from 
0  =  0  to  6  =  2  it  for  the  selected  R.     The  plot  as  recorded  on  the 
Sanborn  Recorder  is  displayed  as  separate  imaginary  (y)  and  real 
(x)  components.     The  polynomial  contour  (P(z)  plane)  which  may 
be  constructed  from  the  x   and  y  components,   is  symmetrical 
about  the  real  axis.       The  real  component  is  a  mirror  image 
about  180     and  the  imaginary  component  crosses  through  zero 
at  180     with  the  180°  to  3  60     portion  being  an  inverted  mirror 
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image.     Therefore,   as  is  indicated  in  Figure  4,   only  the  0     to 
180     portion  is  required. 

0  7T 


(a)    Synthesizer  Component  Traces 


(b)    Polynomial  Contour 
Represented  by 
Synthesizer  Traces 


Figure  4 

The  P(z)  plane  may  be  drawn  from  the  imaginary  and  real 
components,   and  the  number  of  encirclements  of  the  constant 
term  counted.     This,   as  Cauchy's  Residue  Theorem  states,   gives 
the  number  of  roots  that  may  be  found  at  smaller  radii  (R)  than 
that  selected. 

The  major  item  of  interest  is  the  y  (imaginary)  cross-over 
(where  the  y  component  equals  zero)  in  relation  to  the  constant 
term  a     on  the  negative  real  axis.     This  point  is  represented  by 
a  horizontal  line  at  x  =  -  a     on  the  x  (real)  component  plot.     Only 
the  6  at  y  cross-overs  need  be  noted  from  the  imaginary  com- 
ponent plot,   and  from  the  real  component  plot  only  the  fact  of 
whether  the  corresponding  points  are  more  negative  or  less  ne- 
gative than  the  constant  term  need  be  noted  in  order  to  sketch 
rapidly  and  adequately  the  approximate  P(z)  plane.     If  at  a  y 
(imaginary)  cross -over,  the  x  (real)  component  crosses  the  con- 
stant term  line  at  that  point,   roots  occur  at  that  radius,   one  real 
root  if  8  =  ir ,   two   complex  conjugate  roots  if  6  equals  anything  other 
than  n.     This  immediately  gives  the  (ADR)  and  (DR)  of  those  roots, 
since  (ADR)  =  tt-6  and  (DR)  =  cos       (  it  -  Q)3   as  well  as  the  mag- 
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nitude  of  the  roots,   since  R=u     .     The  typical  P  (z)  planes  of 


Figure  5  will  illustrate  this. 


Three  encirclements  of  a 
3    roots    at  u    <  I 


(a)   P(z)  contour  for  R  =  1 


o 


One  encirclement  of  a 

o 

1  root  at  u       <  .  8 

n 

2  complex  roots  at  u   =.8 

r  n 

e  n**t  i=  cos  |=>707 


(b)     P(z)  contour  for  R  =    .8 
Figure  5:    Typical  P(z)  Contours 
The  approximate  P(z)  plane  need  not  be  sketched  at  all 

since  the  encirclement  criteria  is  immediately  apparent  from  the 
x  and  y  component  traces  directly.     If  during  the  interval  between 
two  y  cross-  overs,  the  real  plot  crosses  the  constant  line,   and 
remains  on  the  opposite  side  when  the  y  plot  again  reaches  zero, 
one  encirclement  is  indicated.     That  is,   one  root  on  the  complex 
plane  exists  at  a  radius  less  than  that  for  which  the  plot  was  made, 


This  is  illustrated  in  Figure  6. 
0  * 


■^ 

0                   6 

I         fi 

2        ff 

P(z)  plane 


One  encirclement  between  0  and  6. 
No  encirclements  between  0,  and 
(9?  or  it 

(a)    One  Encirclement  Indicated  by 

Synthesizer  Traces 

Figure  6 


(b)    One  Encirclement  In- 
dicated by  P(z)  Contour 
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It  should  be  noted  that  from  the  component  or   p  (z)  plot  of 
just  a  single  radius,   no  conclusion  should  be  made  as  to  whether 
the  roots  indicated  are  real  or  complex,   nor  estimate  the  (DR) 
of  such  roots  unless  some  additional  information  is  known  about 
the  equation  being  tested  or  unless  the  x  plot  crosses  the  constant 
line  at  or  very  near  the    G   at  which  the  y  cross-over  occurs. 
Caution  should  be  exercised  when  attempting  to  draw  concrete 
conclusions  in  estimating  (DR)  from  plots  of  x  plots  crossing  the 
constant  line  near  a  y  cross -over  until  the   effects  of  R  variation 
for  the  particular  equation  is  fully  appreciated  and  relative  mag- 
nitude of  synthesizer  indications  and  tape  scaling  factors  are  consi- 
dered.    Examples  (I)  and  (2)  demonstrate,   among  other  things, 
how  an  error  in  estimating  (DR)  may  occur  by  considering  only 
the  results  of  a  single  run 

By  making  several  runs  of  constant  R  for  varying  G  on  the 
machine,   a  composite  plot  may  be  constructed  as  is  shown  in 
Example  (3). 

The  magnitude  of  the  roots  may  be  estimated  within  the 
limits  of  the  selected  R's  by  the  criteria  described  above.     Curves 
connecting  points  of  the  various  R  plots  for  equal  angles  (6)  pro- 
duce curves  that  are  actually  constant  6  curves  for  varying  R. 
The  number  of  encirclements  of  the  constant  term  by  these  con  - 
stant  G  curves  give  the  number  of  roots  at  0's  greater  than  the 
value  of  6  for  that  curve.     This  gives  the  limits  on  the  damping 
ratio,   since   t,     -cos      (180°-  0). 

As  is  indicated  in  Figure  7,  the  number  of  encirclements  of 
the  constant  G  curves  is  obtained  by  considering  G  to  be  held  constant 
until  an  infinite  radius  (R)  is  reached  then  permitting  G  to  complete 
its  sweep  to  180°.     The  latter  part  of  the  sweep   (  @  infinite  R)  will 
never  appear  on  any  trace  or  plot  but  must  be  considered  in  order 
to  determine  encirlements, 
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'  R=oo 


R=0 


.0  =  135  0<:R£oo 

135°  £  0<36O°  at  R=oo 


0  =  120 


(a)    Constant  0  curves  for 
Third  Order  Equation 


R  =  Q0 


R=0 


R=» 


0  =  135 


R=0 


0=  120 


(b)    Constant  0  curves  for 
Fourth  Order  Equation 


Figure  7:    Typical  Constant   0-  Variable  R  curves 
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CHAPTER  4 

TECHNIQUES  INVOLVED  IN  THE  TREATMENT  OF  SPECIFIC 

EXAMPLES 


EXAMPLE  I:     Constant  Coefficient  Case 

In  order  to  scale  the  problem  in  a  manner  suitable  for 
setting  into  the  synthesizer,  the  polynomial  in  z  is  written  in 
the  polar  form  as  shown  below. 

The  polar  form  coefficients  are  determined  for  each  of  the 
trial  radii  (R's)  selected.     The  scale  setting  is  then  determined 
so  that  the  largest  coefficient  may  correspond  to  the  largest 
synthesizer  scale,  setting  (1.  0)  in  order  that  the  maximum  syn- 
thesizer accuracy  may  be  obtained.     This  simple  set  df  slide 
rule  ratios  lends  itself  to  tabular  presentation  as  is  indicated 
in  Table  I,   page  38  which  follows  the  examples. 


P  (z)  -  52  =  z5  +  13  z4  +  60  z3  +  120  z2  +  124  z 


=  R5   eJ59  +  l3R4e349+60R3eJ3e 

+  120  R  2    e    J20  +    124    Re    J0 
since      z    =   Re 
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EXAMPLE    I 


Constant  Coefficient  Case 


The  Sanborn  record  of  the  synthesizer  component  traces  are  sketched  and  summarized  below  for  convenience. 
The  values  at  the  important  points  and  cross-overs  are  included  and  the  corresponding    P(z)   plots  sketched  in  order 
that  the  effect  of    R  variation  may  be  readily  observed.     The  conclusions  and  estimates  based  on  the  individual  runs 
are  given  with  the  corresponding  traces  in  order  that  the  technique  used     in  the  investigation  may  be  followed  con- 
veniently. 

P(z)   -  52  =  z' 


0°         47°      92°  127°  148°  w 


5  encirclements;  hence  5  roots  at  R  <  6 
At  least  1  real  root. 

One  complex  pair  may  be  expected  with 
a  d  at 

148°  >  6  >  127° 
or 

127°  >  6  >  92° 


2  roots  @    R  =   V26 

e  -  145° 


c 


35 


(The  run  for   R=6  indicated  the  possibility 
of  roots  (9  6  =  145?  The  other  indication,   i.e. 
127^e>92°did  not  materialize) 

3  encirclements 
3  roots  §  R<  S26 

Due  to  the  position  of  the  loop  in  the  P(z) 
plot  the  only  safe  estimate  is  that  if  there 
is  another  complex  pair,  it  should  occur 
at     100^:  9  <  ir 


0°  51°        100°        145°    w 


R*26 


3  roots    @  R  <  3.  5 

If  there  is  another  complex  pair  at  R<  3.  5 
it  is  more  likely  to  occur  between  111  and  n 
then  between  59°  and  111°.  (The  basis  of  this 
estimate  is  the  position  of  the  constant  term 
and  observing  the  behavior  for  many  other 
equations) 


R=3.5 


93°     130°         ir 


2  roots  @    R  =  vT 

8  --  130°       A;  =  50° 
1  root  (real)  @  R  <  vT 


R=12 

x 
0 
-52 


i 


30 


Example  1    (cont. ) 


1  root  0  R<1.  2 

(The  polynomial  contour  plot  and 
component  traces  indicate  the 
real  root  should  be  very  close  to 
R«l.  2) 


R=I2 


x 
0 
-52 


■52 


I  root  @  R*l 
0=r       A,  *  0 

All  roots  have  been  located 


R-IjO 

x 

0 
-52 


-52 


0° 


No  roots  occur 
at  R  <. 5 

This  checks  with  previous  runs 
since  all  roots  have  been  located 
at  R's  >  .5 


R*5 


The  equation      P(z)  =  z5  +  13 z4  +  60z3  +  120z2  +  124z  +  52 


has  roots  of 


z=  -  5±  j       z«-  l±j  z=  -  I 

,o. 


(  un  =  s/l6,    A;  «  11°)    (un«  VT,A;  =  45°)    (Of'l.A  ,=0) 


The  actual  values  of  the  roots  of  the  equation  agree  very  well  with  those 
indicated  by  the  synthesizer. 


31 


EXAMPLE   2 
Constant  Coefficient  Case 


Example  2  is  of  a  form  very  similar  to   Example  1,  even  to  the  extent  of  an  increasing  magnitude  of  the 
coefficients  with  decreasing  powers  of  the  variable. 


Example    1:      P  (z) 
Example    2:      P  (z) 


z      +  108  z  +  52 


(The  constant  term  (a     =  52  in  both  examples  )  while   greatly  affecting  the  roots  has  no  effect  on 
the  shape  or  magnitude  of  the  component  traces  or  polynomial  contour  plot) 

These  two  examples,   while  similar  in  form,   show  quite  different  behavior  and  indicate  that  an  estimate 
of  B  (or  A      )  based  on  a  single  run  may  be  considerably  in  error. 

P(z)  -  52  =  z5  +  5z4  t  36  z3  +  88  z2  +  108  z 


0°         45°     84°  110°  142°      ir 


5  encirclements, 5  roots  @  R<7 
If  complex  roots  exist,  one  pair 
may  be  expected  @    142>8  >  110 
or  110°>8>84° 


R=7 


2  roots  S  R  =  "^26 
9  =  101°    A      =  79° 
Second  estimate  from  the  R  =  7  run 
is  substantiated.     (Compare  the 
R  =   v26  runs  of  Example  1  and 
Example  2) 

3  roots  S  R  <  "/26 


R=126 


0°  67°  1 36"       j 


3  roots  @  R<3 

If  more  complex  roots  exist,    it 
cannot  be  stated  positively  whether 
they  will  be  at    67?6<l36° 

or     I36°<0<ir 
The  relative  magnitude  of  the  trace 
peaks  indicate  however,  that  6  will 
be  closer  to  136°  than  either  67°  or  » 


2  roots  @  R  =   vT 
6=  135°      A     ■  45° 

I  root  (real)  @  R  <  \/2" 


R^ 
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Example  2    (cont. ) 


-52 


I  root  d 
e  =  ir    A. 


R  =  l 

=  0 


All  5  roots  have  been  located. 
No  indication  of  roots  should  be 
found  at  R's  <l 


T 

-52 


No  roots  @  R<.  5 


R=.5 


The  polonomial: 


P(z)  =  z5  +  5z4  +  36z3  +  88z2  +I08z  +52 


has  roots:  z=-l±5j         z«-l±j  z  =  -l 

(u>n  =  ^/26"   ,    A      =  79°)       (u  n  =   VT,  A,  =45°)     (    a  *  l,A.=0°) 


In  this  example  the  actual  values  of  the  roots  agree  almost  exactly  with 
the  values  indicated  by  the  synthesizer. 
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EXAMPLE  3 


Constant  Coefficient  Case 


P  (z)  -  120  =  z3+l0z2  +  44  z 


-120 


-120 


3  encirclements 
3  roots  @  R<7 

(Possibility  of   one  real  root  is 
indicated  at  R  less  than  but 
close  to  7) 

(  Suspect  one  complex  pair  @ 
85°<e<  140°) 


R=7 


1  root  (real)  @  R  =6 
6  =*      A        =  0° 

2  roots  @  R  <  6 

If  a  complex  pair  exists  it  should  be 
found  in  the  vicinity  of  6  -  139 


R=6 


2  roots  @  R<  5.  25 


R=525 
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Example  3    (cont. ) 


2  roots  @  R  =  4.47 
0  =  116°     A    =  64° 


R=447 


16°  134°       ir 


r 


-120 


No  roots  % 
R  i  3.5 


R=35 


The  polynomial 


P(z)  =  z3  +  10  z2  +  44z  +  120 


has 

roots 

z 

6          z 

=  -  2  ± 

4j 

(wB- 

6' 

•\ 

=  0°)  (w     =  4. 

47, 

The  factored  form 

is  therefore 

P(») 

= 

(z  + 

6)(z  +  2 

+  4j)(z 

+  2 

60°) 


4j) 


The  accuracy  of  9  indicated  by  the  synthesizer  is  very  good  except  when  $  approaches  180 

Figure  8  shows  all  the  polynomial  contours  sketched  above.     It  was  constructed  by  taking 
the  values  from  the  tapes  to  provide  a  series  of  superimposed  polynomial  contours  for  all 
values  of  R. 

Another  series  of  curves  was  constructed  by  connecting  points  of  equal  9  in  order  that 
constant  6  curves  for  varying  R  may  be  obtained.     Cauch's  Residue  Theorem  is  applied  to 
the  constant  6  curves.     The  encirclements  are  counted  and  the  6  values  determine  the  sector 
in  which  roots,  if  any,   are  located. 
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EXAMPLE   3 


3  encirclements 
3  roots  <J» 
e>90° 


/' 


r  =3.5 


Izo" 


r=447 


-0*120° 

encirclement 
I  root  (real) 
0>I2O° 
6  =  ir(reol  root) 


flm[P(z)-l20] 


r  =  5.25* 


P(z)-l20  =  z3+IOz2+44z  =  r3ei3fi+IOr2ei2fl  + 

44reifl 


-10 

-20 

-30 

-40i 

-50 

-60 

-70 

-80 

-90 

•100 

-110 

-120 

-130 

-140 

-150 


10        20       30 
Re[P(z)-l20]  — »- 


r=7 


Figure  9=  Expanded  Portion  of  Region  Near  Origin  in  Figure8  for  Constant  8-  Variable  R. 
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EXAMPLE  4:    Adjustable  Coefficient  Case 

Example  4  illustrates  one  analytical  method  of  investigating 
a  typical  equation  with  adjustable  coefficients  in  order  that  the 
limits  of  the  component  parameters  affecting  the  coefficient  may 
be  determined  which  meet  the  specified  dynamic  conditions.     This 
method  as  explained  in  Example  4  utilizes  the  algebraic  relations 
between  the  adjustable  coefficients  and  the  sums  and  products  of 
the  resultant  roots.     This  method  is  suitable  for  lower4  order 
equations  (second  or  third)  with  one  or  two  adjustable  coefficients. 

EXAMPLE  4 

A  typical  characteristic  equation, of  a  positional  servo- 
mechanism  with  modified  velocity  signal  damping,   is: 


2  7r   r        'p 
v 


2   ^  (ps)  1  +  2  7rr     'p    +     2   V    (ps)    (res) 


'p  +   1    =   0  (I) 


where  'p     =  \L —      ,       P  =      Laplace  Transform  operator 

°n  (PS) 


t  u  =     undamped  natural  frequency 


(Ps) 


of  positional  servomechansim 


t,  .     .  =     damping  ratio  of  positional 

servomechanism 


,                                        t 
r; =    CT  =  — ^ =      non-dimensional  character- 

(ps)  (ps)  istic  time  of  signal  modifier 


39 


Expanding  Equation  (I): 


'P      +     2£  /~o*    +     2  t 


'  (ps)  '  b  (ps)  (res)  '    2  7r  r 


v 


^     /r^o\     /vqo\  O   ir     r 


(ps)  (res)     2  7T  r 


,], 


2  7rr 


v 


(2) 


Or, 


'p3   +(a+c+  b)'p2+   (be  +  I)  "p  +  b  =    0 


(3) 


where 


a  =  2    I 


(ps) 


b  = 


1 


2  7T  r 


c  =    2  C 


(ps) (res) 


For  zero  forced  dynamic  error, 

2  £ 
FDE    =     (ps)  (res)     A  ,    o     for  finite    u 


U) 


n 


(cm) 


n 


(PS) 


where    A 


(cm) 


(  PS) 

(angular  velocity  of  controlled  member)   =£  0 


Thus,    I  ,     w        .      =   c  =  0 
'    b  (ps)  (res) 


and  Equation  (3)  becomes 


'p  3  +  (a  +  b)  'p2  +    'p  +  b  =  0 


(4) 


Substituting    'p 


w 


n 


(PS) 


letting  p    =    (J  e 

n(r) 


J  9 
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and  defining    R  =  L_L       and    z    =    R  e    ^ 

(PS) 

Equation  (4)  is  :    P  (z)  =    z  3  +  (a  +  b)   z  2  +  z  +  b  =  0  (5) 

If  the  roots  of  the  polynomial  in  z    |p  (z)|  are  z ,  ,   z  „  ,   and   z  ~  , 

P(z)  =   (  z-zt  )   (z-z2)  (z-z3)   =  0  (6) 

3  2 

=  z      -(z,  +  Zp+    z  ~  )    z    +(z,    zo+z2    zo  +  zizo)z- 

zl  z2  Z3  =  °  <7> 

For  a  third  order  equation  with  two  complex  roots,  the  roots 
will  have  the  form: 


z ,     =   x,    =  -   a 


z2  =x2+^2  =  •    ^r    un,   .    +    J  un,   ,      ^l"  ^ 

(r)  (r) 


z3    =x3-jy3  =  -  ^r  wn  ,         "   J  w  n ,    .       ^w/ 

(r)  (  r) 


Using  these  relations  and  equating  coefficients  of  like  terms  in 
Equations  (5)  and  (7)  we  have: 


a  +  b  =    a  +  2   £r    u  (8) 

(r) 

I  =   2  £     u  a+  u2  (9) 

(r)  (r) 

b    =    a    J  (10) 

(r) 
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For  the  purpose  of  illustration,   the  specifications  for  the 
dynamic  characteristics  were  taken  to  be 


I)         .4  <  u  <l 

-      n  .   .  - 
(r) 


2) 


t   r  ^'5 


3)     a     and   £     w  >  .  2 

(r) 


Setting    £,        =     .  5  : 


From  Equation  (8)       a  +  b   =   a-  +w 


n 


(r) 


or     a    =  (a  +  b)-u 


n 


(r) 


Using  Equation  (9)  I  -  u 


n 


w 


n 


(r)  (r) 


(a+b) - u 


n 


(r) 


or    u 


n 


(r) 


a+b 


Thus  a    =    a  +  b  - 


a  +  b 


(  a  +  b)      -I 

a+T5 


and     b   =   a*u 


(a  +  b)    -I  I  (a+b)  _    -I 


n 


(r)  a  +  b  (a  +  b)2 


(a+b) 


Choose    a+b,   calculate    b. 


For   b 


db 

max  '     d  '  t 


:■•   + 


(a+b)2         <a  +  b)4 


-(a  +  b)     +  3 
(a+b)4 


=  0 


.-.  (a  +  b)|    b  =     ^^3" 


max 
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and  b  ^aY    =      {sfl)     :  l     =    tt-A-   =    .385 
max  .    iov  3  o\l    3 


Setting  u> 

(r) 


When     ^r=    .5,  a  +  b  =  a  +  ,  4 


a     =  (  a  +  b  )  -  .  4 

l-(.4)2  =  2  (.5)   (.  4  )  |(  a+b)  -.  4 

84 
a+  b  =       ^p    +    .   4    =    2.5 

a  =  2.  5  -  .4     =    2.1 


b  =  (2.1)   (.4)2  =   .336 


When  £r  =  l.  a  +  b  =    a-  +  .  8 


a  =   (a  +  b)  -  .8 

I-    .  42  =   2  (1)  (.4)  [(a+b)-  .  8j 

84 
a  +b   =      ^  +    .  8  =  1.  85 

a   -  1.85  -  .  8    =  1.05 


b   =    (1.05)  (.4)2  =  .168 


Other  points  between  these  limits  are  calculated  similarly  by 
using  other  values  of     .  5  <  t,     <   I 


43 


For   a   =    .  2 


b  =  a  u  =  .  2  w' 


n 


n 


id 


n 


5b 


(r)  (r)  (r) 


1  -  u  =     2  r       u  (.2) 

(r)  (r) 


I  -u 


n 


2  £   *.  u 

^  r      n 


(r) 


.  2 


(r)  C    /I  2  X 

^— L    =    5  (1-u  ) 

n  ,   v  ' 


<r) 


a+  b  =  .  2  +  2  C     w 
^r      n 


=    .2  +  5  (l-u*        ) 
(r)  "(r) 


=    .2  +  5  (l-5b) 
=  .2   +5  -  25b 
a  +  b  =  5.  2  -  25b 


b   = 


5.2-(  a  +b) 


25 


Choose  values  of  a  +  b,   calculate  b. 


For  L      =  I 

'  r 


a+b  =  a-  +   2  cj  ,      l-u  =    2  u 


n 


n 


n, 


(r)  (r)  "(r) 


(a+b)-2u 


n 


(r) 


b  =  «  u 


n 


(r) 


(a+b)  -2  u 


n 


(r) 


(0 


n 


(r) 


l_w 

in      2 

2  u  Wn/r.\ 

n(p)  (r) 
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I  -  u 


n 


b    = 


(r) 


U) 


n 


(r) 


a+b 


1  -u 

■5 &-    +     2u 

2  u  n,   v 

n  (r)  (r) 


I  +  3    u 


n 


(r) 


2  u 


n 


(r) 


Assume    0<  u)         <  I,   calculate  corresponding    a  +  b     and 
r 


b  on  t       =1    line. 

*  r 


From  the  preceding  calculations',   Figure  10  was  constructed, 
which  shows  the  possible  combinations  and  variations  of  the 
parameters  (  a+b)  and  (b)  which  satisfy  the  performance  specifi  - 
cations. 
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The  complex  plane  root    loci    in  Figure  11  may  be  constructed, 
with  the  variables     (a+b)  and  (b),   as  shown  below: 

P  (z)  =  z3  +   (a+b)   z2  +  z    +  b   =    0 

Relating  equations: 

a  +  b    =  a+    2  I    u>  (a) 


r     n 


(r) 


(r)  (r) 

b     =    a   u  (c) 

(r) 


Solving  (b)  for    a-  : 

2 


1  -  U) 

a   =   ill  m 

a  2£       u  *1' 

(r) 


Substituting  in  (a)  : 

2 


l-u 

a  +  b    =-2-TlF    +    2^r%,   ,  <2> 

r      n(r)  <r> 


Solving  for   2  £      u 


J- 


(a+b)  ±  -7  (a+b)    -  4    (l-u„          ) 
2;     y  =  a (3) 


By  assuing     0  <  u  <    1.0,  2^u  can  be  calculated 

n(r)  r    n(r) 
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in  terms  of  a  +  b. 

For  each  assumed  value  of    u  ,  values  of  a  +  b  may  be 

tabulated,   and  the  corresponding    a '  s  and  b  '  s  found  which  satisfy 
the  above  relating  equations. 

These  tables  are  presented  in  Appendix  B  and  the  results 
are  plotted  in  Figure  ll„ 

The   a  =     2  line  was  obtained  by  letting  a  -  .  2  in  Equation  (a), 

Solving  for   L     u  : 

(r) 


I       u  =       <a  +  b)  ~a    =    (a  +  b) 


r       n. 

(r) 


The  points  on  the  a  -  .  2    line  will  lie  at  the  intersection  of 

the   a  +  b  ~  constant  curves  and  the  corresponding  vertical    £      (J 

n(r) 

lines.     This  curve  is  also  plotted  in  Figure  11. 

The  area  enclosed  by  the  a  =  .  2  line,  the   t,      =  .  5  line,   the 

£,      =  1,0  line,   and      4  <  w  is  the  area  in  which  the  values 

r  —     n  .   , 

(r) 

of  a+b   and  b  will  satisfy  the  specifications  of  the  problem. 
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One  might  also  construct  such  a  plot  as  Figure  11   solely 
through  the  use  of  the  synthesizer.     Runs  are  made  at  various 
radii,  (R)  ,   and  for  various  values  of  the  second  coefficient,   (a  +  b), 
and  reading  the  values  of  the  constant  term,   (b),   which  would  give 
roots  at  each  combination  of  (R)   and  (a+b)„ 

This  is  not  recommended,   however,   since  a  very  large  num- 
ber of  runs  would  be  required,   with  the  number  of  runs  required 
being  further  increased  if  there  were  more  adjustable  coefficients, 
Analytical  methods  such  as  this,   or  use  of  the  cubic  charts  in 
Volume  II,  '    '    Instrument  Engineering  would  enable  one  to  construct 
a  coefficient  boundary  plot  in  less  time. 

In  engineering  applications  and  control  system  design  prob- 
lems such  stringent  boundaries  may  not  be  required.     If  the  pri- 
mary requirement  is  selection  of  components  with  parameters  (repre 
sented  by  the  coefficients)  which  will  permit  the  system  perform- 
ance to  be  well  within  the  specifications,   use  of  the  synthesizer 
may  greatly  expedite  selection  of  coefficients  which  appear  to  be 
satisfactory.     Checking  this  selection  on  the  synthesizer  will  insure 
that  the  performance  is  within  specifications. 

In  this  instance  only  two  runs  on  the  synthesizer,   one  at  a 
radius,   (R),  of  about  1,  0  and  one  at  a  radius  of  about  0.  7  for  an 
(a  +  b)  of  around  I.  5  will  be  sufficient  to  give  a  rough  but  adequate 
locus  of  complex  roots  on  the  root  plane.     Two  more  runs  at  an 
(a  +  b)  of  about  l«  7  at  the  same  radii  will  give  another  locus.     These 
four  runs  and  two  resultant  loci  will  be  sufficient  to  give  a  general 
idea  of  the  loci  of  other  selected  values  of  the  coefficient,     Figure 
12  shows  such  an  estimated  plot  which  corresponds  to  the  computed 
plot  of  Figure  11. 
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Approximate  Root  Plane  Plot  of  Example  4 

Figure  12 
From  this  rough  root  plane  locus  it  is  apparent  that  a  value 
of  (a+b)  of  1.  6  with  a  constant  term  (b),  of  0.  2  should  be  one  of 
many  possible  combinations  which  give  a  satisfactory  solution 
well  within  specifications.     There  remains  only  to  set  up  this 
equation  to  check  the  proposed  selections.     Few  trial  and  error 
runs  will  be  required  to  locate  the  roots  accurately,   since  a 
glance  at  the  rough  root  plane  plot  will  show  that  the  complex  pair 
should  be  found  at  an  w     ,  of  about  0.  65  or  0.  7.     A  rudimentary 
knowledge  of  Root  Locus  will  indicate  that  the  real  root  should  be 
found  at  a  radius  of  about  0.  5  since  for  this  equation  the  real  root 
moves  left  from  the  origin  at  twice  the  rate  the  x  component  of 
the  conjugate  pair  moves  right.     The  Root  Locus  method  is  de- 
scribed in  Volume  II,   Instrument  Engineering  *   '. 

As  may  be  seen  from  Figure  10  giving  the  boundary  values 
of  (a+b)  and  (b)  which  yield  acceptable  solutions,  the  three  adjust- 
able parameters,  a,  b,  and  c,   reduce   to  two  adjustable  coefficients, 
It  is  shown  that  boundary  values  of  two  adjustable  coefficients  may 
be  presented  on  a  two-axis  or  two-dimensional  plot.     It  follows  that 
with  three  adjustable  coefficients,  which  may  be  made  up  of  any 
number  of  adjustable  system  parameters,  a  cubic  or  three-dimen- 
sional plot  would  be  required. 


51 


It  is  apparent  that  with  more  than  two  or  three  adjustable 
coefficients,   the  presentation  of  definitive  boundary  conditions 
for  each  coefficient  when  dependent  upon  each  of  the  others   be- 
comes impractical.     It  was  also  found  that  a  method  based  on 
the  algebraic  relation  between  adjustable  coefficients  and  re- 
sultant roots  becomes  unwieldy  for  fourth  order  or  higher 
equations. 

i.e.,      Ist     coefficient    =    I 


n 


z  . 
1 


2  nd    coefficient    =    \ 

i=l 

n 

nth      coefficient    =  .  "L         z  i 

i=0  i 


Accordingly,   it  appears  that  the  most  generally  suitable 
method  of  investigating  the  result  of  adjustable  parameters,   and 
hence  adjustable  coefficients,  will  be  to  set  up  a  series  of  in- 
equalities for  various  conditions  and  observing  the  effect  upon 
the  resultant  performance.     From  these  inequalities  reasonable 
engineering  estimates  may  be  made  for  the  selection  of  various 
sets  of  coefficients,   and  the  resultant  dynamic  characteristics 
that  each  set  should  yield.     These  estimates  are  readily  checked 
on  the  machine  in  the  manner  given  for  the  fixed  coefficient  case 
to  insure  that  they  do  in  fact  give  performances  close  to  that  which 
was  anticipated.     These  checks  performed  on  the  machine  give, 
in  addition  to  rather  precise  dynamic  characteristics  for  the 
coefficients  selected,   a  good  indication  of  the  effect  of  certain 
other  variations  in  selected  coefficients.     This  technique  will  be 
illustrated  in  Example  5. 
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CHAPTER  5 
MODIFICATIONS  INSTALLED  ON  THE  FOURIER  SYNTHESIZER 


The  considerable  labor  required  in  cross  plotting  from  the 
x  and  y  tapes  to  the  P  (z)  plane  may  be  eliminated  by  the  use  of  a 
complex  -  vector  adder  connected  to  x  and  y  signal  leads.     A 
servo  driven  two-axis  recorder  was  used  for  this  purpose.    Timed 
pulses  were   introduced    on  the  x  and/or  y  axis  to  provide  an  in- 
dication of  selected  values  of  8.     The  harmonic  synthesizer  sweeps 
out  from  0°  to  360     at  a  uniform  rate  and  6  is  therefore  propor- 
tional to  time. 

A  more  satisfactory  method  of  observing  and  recording  the 
polynomial  contour  (P  (z)  plane)  plot  was  found  by  utilizing  a 
standard  cathode  ray  oscilloscope  in  place  of  the  two-axis  recorder, 
and  a  Polaroid  Camera  for  recording.     The  associated  circuitry 
and  additional  components  required  to  modify  the  harmonic  syn- 
thesizer to  provide  a  suitable  computer  are  described  in  Appen- 
dix A.     This  resulted  in  a  much  smoother  and  more  accurate  plot 
since  the  mechanical  and  servo  response  difficulties  inherent  in 
the  two  axis  recorder  were  eliminated.     Also,   since  the  area  of 
primary  interest  on  the  P  (z)  plane  is  in  the  vicinity  of  the  x  axis 
cross-over  point  (where  y  =  0),  this  area  can  be  enlarged  to  cover 
the  entire  scope .     Great  latitude  is  possible  in  scaling  both  the 
synthesizer  and  scope,   permitting  substantial  flexibility  in  select- 
ing the  regions  to  be  more  closely  examined.     Through  the  use  of 
the  Polariod  camera  and  multi-exposure,   as  many  runs  as  neces- 
sary can  be  superimposed  on  one  record  for  comparison.     The 
sweep  speed  of  the  trace  is  slow  enough  so  that  it  is  practicable  to 
follow  the  spot  with  a  grease  pencil  directly  on  the  face  of  the 
scope  or  on  tracing  paper.     Thus  the  camera  may  not  be  required 
unless  a  permanent  accurate  record  is  required. 
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This  superposition  of  record  technique  is  the  only  prac- 
tical method  by  which  curves  for  constant   9   and  varying  R  can 
be  obtained,   since  the  synthesizer  is  inherently  a  constant  R 
and  variable   9  device. 

The  set  of  constant    9   curves  are  obtained  merely  by  con  - 
necting  points  of  equal   9   on  the  several  constant  R  curves.     For 
these  records  the  synthesizer  is  not  scaled  for  maximum  amplitude 
for  each  run  of  different  radii,   but  scaled  as  shown  in  Example 
5  to  provide  the  actual  relationship  of  each  radius  plot  to  all  the 
others. 

6  indication  was  provided  by  fabrication  and  installation 
of  a  cam  on  the  synthesizer  drive  shaft  as  described  in  Appendix 
A.     A  microswitch  actuated  by  the  cam  follower  was  installed  on 
the  synthesizer  and  the  signal  produced  was  used  to  blank  out  the 
scope  for  about  4     ,  thus  providing  a   9  indication.     Cam  teeth 
were  cut  at   6   =   90°,    120°,    135°,   and  150°,   and  the  cam  was  cut 
to  blank  out  the  scope  from  180°  to  3  60°,   since  the  trace  is 
symmetrical  about  the  x  axis  and  the  last  180°  would  contribute 
nothing  further. 

The  x  and  y  axes  on  the  scope  trace  may  be  obtained  simply 
by  turning  off  the  y  and  x  scope  amplifiers  respectively.     The 
x  and  y  axes,   being  linear  scales,    may  be  calibrated  in  a  manner 
identical  to  the  method  of  calibrating  the  component  tapes  or  as 
shown  in  Example  5.     Thus,  the  x  axis  location  of  the  constant 
term,   (encirclement  point)  may  be  determined.     The  dynamic 
characteristics  or  roots  of  an  equation  with  constant  coefficients 
may  be  found  to  a  satisfactory  degree  of  accuracy  quickly  and 
easily  without  using  large  amounts  of  Sanborn  Recorder  tape. 
Various  radii  (R)  runs  are  scaled  and  set  into  the  synthesizer. 
Observing  each  trace  for  encirclements  and  the  x  axis  cross  - 
overs  relative  to  the  constant  term  on  the  x  axis  serves  as  a  good 
estimate  for  the  radius  of  the  next  run.     Scaling  and  setting  values 
into  the  machine  requires  very  little  time  or  effort.     All  roots  of 
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the  equation  may  be  found  in  this  manner.     A  root  is  indicated 
when  the  trace  passes  through  (or  near,  depending  on  the  degree 
of  accuracy  desired)  the  encirclement  point.     The  radius  (R) 
used  gives  the  u      of  the  root,  and  the  6    marks  indicate  the 
polar  angle  of  the  root,   measured  counter-clockwise  from    n 
on  the  root  plane     [(ADR)   =  ir  -  0j  ,  as  well  as  the  damping  ratio 

\t,       =     cos    (ADR)j    .     No  Sanborn    tape  need  be  used  unless 
desired.     When  a  more  precise  value  of  (ADR)  is  required  than 
that  obtainable    from  the  five   8  indications  provided  by  the  cam, 
a  single  tape  can  be  run  for  the  root  value  of  the  radius  R  and 
a  more  accurate  value  of  6  obtained. 

For  higher  order  equations,  fifth  or  higher,   it  may  be 
desirable  to  enlarge  the  scope  presentation,  particularly  when 
searching  for  roots  at  the  smaller  radii,  so  that  the  immediate 
vicinity  of  the  encirclement  point  occupies  the  entire  presentation  . 
When  the  origin  no  longer  appears  at  the  far  right  side  of  the  scope, 
as  with  an  enlarged  presentation,   it  may  be  necessary  to  re  -  run 
previously  scaled  plots  where  the  x   cross  overs  are  known  in 
order  to  re- scale  the  portion  of  the  negative  real  axis  that  ap- 
pears on  the  presentation.     With  some  higher  order  equations  it 
may  be  advantageous  to  divide  through  by  the  constant  term  prior 
to  making  any  runs,  thus  shifting  the  origin  to  the  constant  term 
on  the  negative  real  axis.     Thus  the  origin  becomes  the  encircle- 
ment point  and  calibrating  the  x  axis  is  no  longer  a  problem. 

The  principles,   procedures  and  equipment  just  discussed 
are  particularly  well  adaptable  to  equations  with  fixed  coefficients 
and  an  adjustable  constant  term.     Problems  of  this  type  arise 
frequently,  particularly  in  control  system  and  servo  problems 
involving  an  adjustable  open  loop  sensitivity.     Runs  are  made  as 
previously  discussed  for  a  few  radii  and  the  w      and  (ADR)  at 
x  cross -overs  are  plotted  on  the  root  plane.     The  value  of  the 
x   axis  intercept     is  the  constant  term  for  the  equations  thus  formed 
with  roots  at  the  u        and    t,   indicated  by  the  cross- over.      Or- 
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dmarily  just  a  few  such  runs,   with  their  corresponding  roots  and 
value  of  the  constant  term  will  be  sufficient  to  fair  in  a  reason- 
ably accurate  root  locus  on  the  root  plane  plot.     The  effect  on  the 
magnitude  and  (ADR)  of  the  roots  of  varying  the  constant  term 
(or  open  loop  gain)  may  be  determined  from  the  tape  or  photo 
record. 

The  analysis  of  the  case  where  all  the  coefficients  and 
the  constant  term  are  adjustable  will  be  illustrated  by  an  example 
in  the  next  chapter,     This  case  arises  when  system  parameter 
specifications  are  to  be  set  which  will  enable  the  resultant  system 
to  meet  the  required  dynamic  characteristics. 
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CHAPTER  6 
ILLUSTRATION  OF  THE  ANALYSIS  TECHNIQUE 

EXAMPLE  5:    Adjustable  Coefficient  Case 

A  typical  performance  equation  of  a  lead  modified  positional 
servomechanism  is 


»p'+«dPt£*?» 


"! 


(PS) 


Acm  ■  irfr  Win  +  Vtf)  (l' 


where     I     =     equivalent  moment  of  inertia,  referred  to  the 
controlled  member 

c,     =     damping  coefficient,  referred  to  the  controlled 
member 

1+  y  t  p     _     performance  function  of  the  lead  modifier  net- 
1  +  t  p 

work,  where  7  <  I. 

t      =     characteristic  time  of  the  lead  modifier 

S,      »     =      sensitivity  of  the  positional  servomechanism 
for  angle  in,  torque  out. 


M 


.   ,»     =    interference  torque 


A  =     angular  position  of  the  controlled  member 

cm  &  r 


A.  =    command  input  angle 


p  =   j7  =    Laplace  Transform  operator 
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If  we  define: 


I 


2  £ 


7  S  (PS)  J 


(Mi 


7  S 


u 


(PS) 


(PS)  n(ps) 


LJ 


and     r 


n 


(PS) 


n 


(PS) 


Equation  (I)  becomes 


r, 


^2 


P     +  2  C 


,  H-T27rr'p 

(ps)       p        7  (1+2  tt  r  '  p)    rt  out 


1+  7  2  7T  r  '  p    .  f9 

7(1+2  7T  r  'p)        in       *    ' 


so        'p      + 


2t/no»   + 


1 


(ps)        2irr 


'p2  + 


2£ 


(PS) 


27rr 


+    I 


r       7     2?rr  \      out 


7 


2?rr 


+    7  'p 


A. 


in 


(3) 


The  characteristic  equation  is: 


2C 


P3  +  (2  5/^+9^")    'P2    +(-M  +l)'P+4-^r  =  0  (4) 


(ps)      27rr 


2?rr 


7    2?rr 


i0  I 

let        'p  =  z  =  R   e J     ,     2  L ,      ,     ^    a,     x —      =  b,    7    =    c 

r  '  (PS)  25TT  ' 


Then    P  (z)  -    -  =    z3  +    (a+b)   z2  +  (ab+l)  z 


(5) 


or  P(R,0)  »  —  =  r3    eJ30    +    (a+b)   R2   e  J26  +  (ab+l)   R  e  J6      (6) 
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Taking  reasonably  broad  limitations  based  on  servo  principles  as 
a  starting  point  for  the  analysis: 

.5'^;,        <l.  5  .  05^r<.5  ,1*7*  I 

Hence:    I  <  a  <  3  .  3  <i  b  ^  3  .  I  *  c  *  I 

and        .  3*  ab  ^  9 

Therefore:    I.  3<a+b^6  I.  3<ab  +  l^l0  .   I*  c  <  I 

Values  of  the  second  coefficient,  a  +  b,  were  picked  which 
covered  the  range  of  variation  of  a  +  b;  i.  e.        a+b  =  I.  3,   3.5,   and  6. 

However,  choosing  a  value  for  a  +  b  restricts  the  possible 
variation  of  the  third  coefficient,  ab+l,  since  b  must  remain  a 
real  number  physically. 

Thus,   if  a+b  =  I.  3, 

a  =  I.  3  -b 

and    ab+l  =  x    =    (I.  3-b)  b  +  l 

=  -b2  +  I.  3  b  +  I 

Transposing  and  solving  for  b  , 


b        l.3±    Vl769  -4(x-l)  (7) 

2 


The  maximum  value  for  x  for  which  b  remains  a  real  number 
is  obtained  by  setting  the  quantity  under  the  radical  equal  to  zero 
and  solving  for  x; 

4  (x  -  I)  =  I.  69 

x   =169  it 
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Therefore,  1.3   <  ab+l  <  1.42         when    a+b  =  l.  3 

The  value  of  b  is  determined  from  Equation  (7)  by  the 
substitution    x  =  1.42. 

Similarly,   the  permissible  variations  of  ab+l  when  a+b  =  3.  5 
and  a+b=6  were  determined  and  found  to  be 

I.  3  <  ab+l<  4.  06  when  a+b  =  3.5 

and     1 .  3  <  ab+  I  <  10  when  a+b  =  6 

On  the  basis  of  the  preceding  estimates,   runs  were  made 
for  various  R's  with  the  following  adjustable  coefficient  settings: 


Run 

a+b 

ab+l 

_b 
c 

l-l 

1.3 

1.3 

.  3<-<  3 

~~  c  — 

1-2' 

1.  3 

1.42 

.  65<-<  6.5 
c 

2-1 

3.5 

1.3 

.3   <-  <  .  9 

-  c  - 

2-2' 

3.5 

4.  06 

1.75<~<  17.5 

—   C  "" 

3-1 

6 

1.3 

.3<-<  .5 
c 

3-2 

6 

10 

3<-  <  30 
c 

These  runs  were  photographed  and  are  presented  on  the 
succeeding  pages.     The  conclusions  drawn  from  the  examination 
of  the  photographs  appear  with  the  records. 
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Figures    13,   14,   15,   16,   17  and  18  For  Example  5   Follow. 
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For  any  general  case  where  it  is  desired  to  select  coeffi  - 
cients  specifying  component  parameters,    the  effect  on  the  dynamic 
characteristics  of  varying  any  coefficient  may  be  readily  observed. 
This  is  shown  in  the  following  series  of  photographs.    Figs.   19 

and  20. 

The  equation  used  is  of  the  form: 

P(z)    =    z3  +  Az2  +  Bz  +  C    =    0 

A  was  taken  as  1,  3,  3.  5  and  6„ 
B  was  taken  as  1.  3,  10  and  1  9, 
C  is  unrestricted  and  to  be  selected 

from  whatever  specifications  are 

imposed. 

It  will  be  observed  that  this  is  simply  an  extension  of  the 
estimated  boundaries  of  coefficients  of  the  previous  specific 
problem,     Some  of  the  same  runs  were  used  to  demonstrate  that 
there  is  nothing  unique  about  any  one  set  of  runs.     They  might 
be  equally  applicable  to  any  one  of  a  variety  of  problems. 

For  the  purposes  of  illustration,   one  set  of  runs  will  be 
analyzed  for  possible  conclusions  relating  the  dynamic  charac- 
teristics and  the  coefficients.     All  other  sets  for  any  problem 
are  handled  in  a  similar  manner. 

Photograph  3-3  representing 

P(z)    =    z3  +  6z2  +19z   +  C    =    0 

is  sketched  in  Fig.  21  for  clarity  although  constant    6  curves  may 
be  drawn  directly  on  the  photograph  if  desired 


68 


Figures  19  and  20  which  follow  show  the  effect  of  the 
variation  of  the  equation  coefficients. 
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P(z)  =  z3+Az2  +  Bz-i-C 


Analysis  from 
Constant  R  curves 


10^*15    lcomplex  pair  5  <w<6  \ 


1  real  root 


5*<*?1.  0 


\ 


\ 


15?C?20    1  complex  pair    4?w    ?5 


1  real  root 


Ka<3 


25*C?65    1  complex  pair  3<wn?4 

1  real  root  3  ?o<5 

65?C?140   lcomplex  pair  4?w    ?5 
r  n 

1  real  root  65<or<4 


\    e=i2oe 

\ADR=90 

y.o 

\ 
1 


Analysis  from 
Constant  8  curves 

0<C<20      1  complex  pair  30°  <ADR< 45° 
1  real  root  ADR  =    0 

C  =   20  1  complex  pair  ADR  =  30°(\vn  =  4) 

1  real  root  (a  -  1.5) 

20"?C*30    1  complex  pair  3<P<ADR<45° 
1  real  root 

30?C*40    1  complex  pair  45°  <  ADR  <  60° 
1  real  root 

40*C?250  lcomplex  pair  60°<ADR<90° 
1  real  root 

250  <  C  1  complex  pair  ADR  >  90° 

1  real  root 

Unstable  since  one  complex  pair 
occurs  in  right  half  plane 


Fig.  21      Polynomial  Contours  from  Photographs  3-3  Constant  R  and  Constant  6  Curves. 
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CHAPTER  7 
CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS: 

In  this  thesis,  the  mapping  theorem  was  used  in  estimating 
the  dynamic  characteristics  of  physical  systems  from  the  char- 
acteristic equation  of  the  performance  function.     Modification  of 
a  Fourier  Harmonic  Synthesizer  greatly  facilitated  the  application 
of  the  theorem,   and  enabled  the  polynomial  contour  to  be  obtained 
directly.     The  presentation  of  the  polynomial  contour  on  a  stand- 
ard cathode  ray  oscilloscope  made  it  possible  to  apply  the  encir- 
clement criteria  conveniently  and  with  considerable  accuracy. 

It  was  found  that  input  signal  levels  too  low  to  be  indicated 
by  the  Sanborn  Recorder  pens  still  produced  accurate  polynomial 
contours  on  the  scope  presentation. 

The  difference  between  the  computed  polynomial  contour 
check  points  and  the  corresponding  points  on  the  photo  records 
was  not  measurable. 

By  expanding  the  scope  presentation,   any  portion  of  the 
polynomial  contour  may  be  closely  examined.     This  procedure  is 
useful  when  investigating  critical,   or  limiting  values  of  the  charac 
teristic  equation  coefficients. 

Examples  have  been  included  which  indicate  that  the  tech- 
niques presented  are  applicable  either  to  polynomials  with  fixed 
coefficients  or  to  polynomials  whose  coefficients  are  adjustable. 
In  the  case  of  adjustable  coefficients,  the  method  presented  pro- 
vides for  the  possibility  of  setting  the  limits  on  the  coefficients 
which  will  allow  the  roots  of  the  system  characteristic  equation 
to  remain  within  a  specified  sector  on  the  root  complex  plane. 
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This  will  allow  the  dynamic  characteristics  of  the  system  to  meet 
the  prescribed  specifications.     The  performance  functions  of  the 
system  components  are  also  determined  and  are  represented  by 
by  the  coefficients.     This  knowledge  would  be  very  helpful  to  the 
preliminary  designer. 

Although  the  examples  presented  here  are  primarily  chosen 
from  servo  problems,  the  technique  is  equally  applicable  to  air- 
craft stability  problems  or  to  other  problems  where  the  dynamics 
determined  by  the  characteristic  equation  are  of  interest. 

RECOMMENDATIONS: 

For  the  purpose  of  illustrating  the  analysis  techniques  for 
the  adjustable  coefficient  case,   a  typical  problem  involving  an 
equation  of  the  third  degree  was  chosen  (Example  5).     It  is  recom- 
mended that  these  procedures  be  applied  to  actual  physical  prob- 
lems of  fourth  degree  or  higher.     The  results  of  this  investigation 
indicate  that  the  techniques  developed  apply  equally  well  to  equa- 
tions of  higher  order  without  excessively  increasing  the  amount 
of  computation  required. 

It  is  further  recommended  that  the  feasibility  of  inverting 
this  technique  be  studied.     The  shape  of  the  polynomial  contour 
is  readily  adjustable  by  the  synthesizer  potentiometers  while  the 
machine  is  in  operation.     Once  the  relation  between  the  polynomial 
plane  contour  and  the  root  plane  is  fully  appreciated,   it  should  be 
possible  to  adjust  the  synthesizer  to  produce  polynomial  contours 
representing  the  desired  dynamic  characteristics.     By  starting 
with  the  potentiometer  settings,   and  inverting  the  scaling  tech- 
nique, the  equation  coefficients  providing  the  specified  dynamic 
characteristics  may  be  determined. 

If  such  a  technique  can  be  developed,  it  should  facilitate 
further  the  problem  of  estimation  of  the  dynamic  characteristics 
of  physical  systems. 
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APPENDIX  A 

IMPLEMENTATION  OF  THE  HARMONIC  SYNTHESIZER 

MODIFICATIONS 


The  harmonic  synthesizer  used  in  obtaining  data  for  this 
thesis  was  constructed  by  the  Instrumentation  Laboratory  of  the 

Massachusetts  Institute  of  Technology  .     The  complete  circuit 

(5) 
diagrams  of  the  synthesizer  are  available  at  the  Laboratory.   'The 

mechanics  of  the  synthesizer  operation  will  be  outlined  in  order  to 

clarify  the  reasons  for  the  modifications  made. 

This  synthesizer  is  composed  of  the  fundamental  and  23 
harmonic  components.     Figure  A  -  1  illustrates  the  operation  of 
the  fundamental  component  of  the  synthesizer.     An  electric  motor, 
through  a  gear  train,   furnishes  the  drive  for  the  rotation  of  the 
rotor  coils,   which  are  90°  apart  electrically.     The  moving  coils 
are  excited  by  a  400  ~  voltage  in  the  stator  coils.     The  signal  vol- 
tages are  picked  off  at  the  summing  resisters,   amplified,   and  sent 
to  a  Sanborn  two-channel  recorder.      Since  the  rotors  for  each 
component  are    90     apart  electrically,  the  resulting  signal  voltages 
will  have  a  phase  difference  of  90°  and  thus  correspond  to  cos    8 
and  sin  Qt  respectively,   as  the  motor  turns  the    rotor  shaft  through 
one  complete  revolution  of  360°.     Multiplication  by  R  is  done  by 
the  potentiometers.     Thus,  the  output  voltages  of  the  fundamental 
component  in  the  synthesizer  correspond  to  R,    cos  6  and  R.   sin  6, 
or  to  the  real  and  imaginary  terms,   respectively,   of 

i  R 
z    =   R,    e    J      =    R.    cos    8  +  j   R,    sin  G 


Similarly,   if  the  rotor  coils  of  the  2      harmonic  component 
have  twice  the  angular  velocity  of  the  fundamental,  the  output 
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voltages  of  this  component  correspond  to 

R9    e    J  26    =  R     cos   20  +   j    RQ    sin    2  e 


z 


and  Z3  =    R3  e    ^  ^    =  R3  cos    3  6  +  J  R3  sin  30'   etc' 


2  3 

If  R,  =  BR,  R2    =  AR     and    Rg     =   R    ,  then  the  sum  of  the 

first  three  components  of  the  synthesizer  is 


P  ( R,  0 )  =  (R3  cos  3  0  +  AR2  cos  2  0  +  B  R  cos  6  ) 


3  2 

+  j    (  R    sin  3  6  +  AR      sin  2  6  +  B  R  sin  0  ) 


Thus,  the  real  and  imaginary  parts  of  any  complex  harmonic 
function  (within  the  limits  of  the  machine )  may  be  readily  obtained 
by  expressing  it  as  a  polynominal  in  polar  form,  such  as 

P(R,0)-a     =Rne  Jn0  +  a      .   R  n  ~l  e  *  <n"  1  >  6     +  .    .    .    . 
x        '       o  n- 1 

....   +  a2  R2  e    J20  +  &l  Re  ^ 

The  characteristic  equation  of  any  physical  system  may  be 
expressed  in  this  manner. 


The  real  and  imaginary  parts  of  the  function  -  henceforth 
referred  to  as  the  x   and  y  signals,   respectively  -  were  applied 
to  the  x  and  y   inputs  of  a  standard  Dumont  cathode  ray  oscilloscope 
which  was  fitted  with  a  Polaroid    Land  Camera 

A  z  -   dot  generator  was  designed  in  order  to  provide  scope 
illumination  at  the  moments  when  the  400  "  excitation  signal  was 
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at  its  maximum,   as  shown  in  Figure  A  -  2. 


Figure  A-2:     Signal  Voltage 
with  Z-Dot  Generation 


These  pulses 
appear  on  the  face 
of  the  scope  at  the 
rate  of  400/  second 
and  appear  to  be  a 
a  steady  dot,  which 
traces  the  complex- 
vector  sum  of  the 
x   and  y  signal  vol- 
tages.    Thus,   a 
visual  represen- 
tation of  the  poly- 
nominal  contour  is 
presented,  which  is 

the  first  requirement  for  the  application  of  the  encirclement  cri  - 

terion. 

Finally,   a  cam -actuated  microswitch  was  installed  on  the 
rotor  shaft  of  the  fundamental  component,  which  provided  timing 
indications  corresponding  to  certain  values  of  6..     The  final   180 
of  the  dot  travel,   being  symmetrical  to  the  first  180°  ,  was  blanked 
out.     Figure  A  -  3  shows  the  Fourier  Harmonic  Synthesizer  after 
modification,   and  Figures  A-4,5,   and  6  are  the  circuit  diagrams 
of  the  x  and  y  signal  amplifiers  and  the  z-dot  generator.     Figure  7 
is  a  sketch  of  the  cam. 

In  operation,   as  the  dot  traced  out  the  plot  of  the  character- 
istic equation  polynomial  ,  the  microswitch  cam  follower  opened 
contacts,  thus  blanking  the  scope,  duringthe  indicated  intervals. 
The  scops  camera,   set  for  time  exposure,   photographed  the  trace 
for  a  permanent  record. 
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APPENDIX  B 


Tabulation  of  Numerical  Values  Used  in  the  Determination 
of  the  Effect  of  Varying  the  Coefficients  in  the  Characteristic 
Equation  of  Example  4: 


P(z)    =  z3  +  (a+b)  z2  +  z+  b 


=  z3  +  (o+2Crwr)   z2  +  (2  i  r  ur  a  +  u  2  )  z  ■+    or  co   2  -0 


where       a+b    -    or    +  2  £        u 

^  r      r 


I     =    2  t  u   o+    w 
^r  r  r 


b     =  a     (j 
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APPENDIX  C 


Tabulation  of  Fourier  Harmonic  Synthesizer  Scale  Settings  Used 
in  Example  5. 


The  basic  equation  used  in  this  example  is 

P(z)  =  z3  +  (a+b)  z2  +  (ab+l)  z  +  — 

i  0 
or,  with     z  =  Re  J 

P(R,0)--|    =    R3  e  j30  +    (a+b)  R2  e    j20    +  (ab+l)  Re  i6 

At    0  =0,   e  J0  =    e  J°  =  I  +  jO 

i20         i0       . 

eJ     =    e  J     =    I  +  j  0 

j30  i0 

eJ      =    e  J     =  I  +  jO 


and  the  real  part  of  P  (R,  0)  -- 


F    =    R3  +  (a+b)  R2  +  (ab+l)  R 
0  =  0    ° 


Similarly,   at    0  =  it,     e  ^     =    e  J  *  -  -  1  +  jO 

eJ       =eJ  1  +  jO 


J30  i3vr 

e  J       =  e  J        =  -  i  +  jo 


and  the  real  part  of  P  (R,  0)  -  |  =  F^  =  -R3+  (a+b)  R2-(ab+l)  R 


0  =it 


Since,   in  this  example,   R  =6,  the  largest  potentiometer 
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3 
setting,   So  will  correspond  to  R  throughout.     That  is, 

max, 

3  3 

division  of  all  numerical  coefficients  by  R  =  6     =  216  gives 

J      max  s 

the  scale  setting  for  the  potentiometer  in  each  run. 

These  tabulations  are  carried  out  in  the  following  tables, 
as  indicated. 


Sample  calculation: 

R  =  6     R2  =  36  R3    =    216 

1.3R2    =    1.3  (36)  =  46.8 

1.3R      *   1.3  (6)    -  7.8 

S3  =  2\6  r3    =  216"  *2l6)  =  L  00°  Fo  =  2l6+  46,  8  +  7l  8  =  2T0-  6 

S2  =  2l6  (13r2)  =  Tle^46,8^0'217       F7r  =  "216  +  46.8  -  7.8  =  -177 
Sl   =216(l'3R)  =  IT6    <7-8)  =  00361 
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TABLE   C  -  I 


Table 

Photo  No. 

Scope  Settings 
(Gain  -  I) 
Y  Amp.            X  Amp. 

Multiply 
Sx    S2  and  Sg  by 

C-2 

l-la 

29                        61.  5 

C-2 

l-lb 

29                        61.5 

100 

C-3 

1-2'a 

24                        55 

C-3 

l-2'b 

24                        55 

20 

C-4 

1=2  a 

18                         44 

C-4 

1-2  b 

18                         44 

10 

C-5 

1-3  a 

22                        48 

C-5 

I -3b 

22                        48 

5 

C-6 

2-la 

30                        64 

C-6 

2-lb 

30                        64 

10 

C-7 

2-2'a 

24                        55 

C-7 

2-2'b 

21                        39 

10 

C-8 

2-2a 

29                        61 

C-8 

2 -2b 

29                        61 

10 

C-9 

2-3a 

30                        64 

C-9 

2-3  b 

20                        45 

20 

C-10 

3-1  a 

20                       37 

C-10 

3-lb 

20                        37 

5 

C-U 

3-2a 

30                        59 

C-U 

3-2b 

30                       59 

20 

C-12 

3-3a 

30                       61 

C-12 

3-3b 

20                       40 

10 
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